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DYNAMICAL FRICTION FORCE EXERTED ON SPHERICAL BODIES 
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RESUMEN 

Siguiendo un enfoque mecanico-ondular calculamos la fuerza de arrastre ejercida por un sistema homogeneo 
e infinito de estrellas de fondo sobre pertubador mientras este se mueve a traves del sistema. Recuperamos 
la formula clasica para la fuerza de friction (FF) derivada por Chandrasekhar, pero con un logaritmo de 
Coulomb modificado. Al estimar la FF ejercida sobre una esfera de Plummer y un perturbador que posee un 
pcrfil tipo Hcrnquist, consideramos un intervalo de modelos que abarca toda distribution plausible de satelitcs 
galacticos. Se muestra que la configuration del perturbador afecta unicamente la forma exacta del logaritmo 
de Coulomb. Tal logaritmo converge a pequenas escalas porque los encucntros entrc la partfcula dc prucba y 
las estrellas de fondo cuyos parametros de impacto son inferiorcs al tamano del perturbador masivo rcsultan 
inefiticntes. Comprobamos asf los resultados previos basados en la aproximacion dc impulso dc pequenas 
deflecciones angulares. 

ABSTRACT 

Following a wave-mechanical treatment we calculate the drag force exerted by an infinite homogeneous back- 
ground of stars on a perturber as this makes its way through the system. We recover Chandrasekhar's classical 
dynamical friction (DF) law with a modified Coulomb logarithm. We take into account a range of models 
that encompasses all plausible density distributions for satellite galaxies by considering the DF exerted on a 
Plummer sphere and a perturber having a Hcrnquist profile. It is shown that the shape of the perturber affects 
only the exact form of the Coulomb logarithm. The latter converges on small scales, because encounters of the 
test and field stars with impact parameters less than the size of the massive perturber become inefficient . We 
confirm this way earlier results based on the impulse approximation of small angle scatterings. 
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1. INTRODUCTION 

The process of dynamical friction (DF) is one 
of the most classical and fundamental problems en- 
countered in the description of the evolution of al- 
most all astrophysical systems. From the critical mo- 
mentum exchange in a protoplanet-protoplanetary 
disk set up, passing through the problem of satel- 
lites in galaxies to galaxies in large clusters, proper 
understanding of DF is a prerequisite to more am- 
bitious attempts at constructing physically justified 
models. 

In his seminal paper Chandrasekharl ( 19431 ) en- 
visaged the scenario of a sequence of consecutive 
gravitational two-body encounters of test a nd field 
stars in order to calculate the drag force (cf. Irlenon 
19731 for a modern presentation). In particular, ap- 
plication of DF to calculate the rate of a sinking 
satellite has received special attention, and efforts 



have been made to include more general background 
distrib u tions . However, as for the perturber itself, 
White! (|l976h has been the only one to consider a 
more realistic finite-size perturber to compute ana- 
lytically the DF based on an impulse-approximation 
approach. His main result was a modification of the 
Coulomb logarithm so that it does not diverge any- 
more at small scales, because gravitational encoun- 
ters at impact parameters smaller than the size of 
the perturbing body become ineffective. Here we 
rigorously calculate the drag force exerte d on differ- 
ent bo dies following the ap proach of both lMarochnik 
(|l968h and lKalnaisf (|l972T ) who determined the "po- 
larization cloud" created in the background medium 
as a massive object was making its way through the 
system. This method can be simply understood as 
linear a nd angular momentum exchange in stellar 
systems (iLvnden-Bell fc Kahiaidll972l |Pekkeij|l974 , 



As- 



1 Astronomischcs Rcchcn-Institut am Zentrum fiir 
tronomie der Universitat Heidelberg, 
Monchhofstrafic 12-14, 69120 Heidelberg, Germany 

2 Fellow of International Max-Planck Research School for 
Astronomy and Cosmic Physics, Heidelberg (csquivel@ari.uni- 
heidelberg.de). 



Tremaine fc Weinberg! Il984l , iFuchsl 12004) , and has 
been extensively used in plasma physics (cf. IStix 
1992). It is shown in section 4 below that the shape 



of the perturber affects only the exact form of the 
Coulomb logarithm. As concrete examples we cal- 
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culate the drag force exerted on a Plummer sphere 
and on a s phere with the density distribution of a 
Hernquistl (jl990t ) profile, respectively, and compare 
this with the drag force exerted on a point mass. 

2. A WAVE-MECHANICAL TREATMENT 

We assume an infinite homogenous distribution 
of field stars on isotropic straight-line orbits. The 
response of the system of background stars to the 
perturbation due to a massive perturber is deter- 
mined by solving the linearized Boltzmann equation 



dfi 
dt 



3 
i=l 



dfx 3$i df 



cbq 9xi <9vj 



= 



(1) 



where $i denotes the gravitational potential of the 
perturber and f± and /o are the perturbed and un- 
perturbed distribution function of the field stars in 
phase space. By Fourier-transforming both the per- 
turbations /i and the potential <&i (whose form are 
fu),k'j ^w.k expi[cjt + k • x] ; where u> and k denote 
the frequency and wave vector of the Fourier compo- 
nents) the solution of Boltzmann equation is greatly 
facilitated. Without loss of generality the spatial 
coordinates Xj and the corresponding velocity com- 
ponents v; can be oriented with one axis parallel to 
the direction of the wave vector. The Boltzmann 
equation |T]) takes then the form 



dv 







(2) 



with k = |k| and v denoting the velocity component 
parallel to k. Equation © has been integrated over 
the two velocity components perpendicular to k. In 
the following we assume for the field stars always a 
Gaussian velocity distribution function, to find the 
solution 



kv 



w + kv ^/2tt(j 



(3) 



where denotes the spatial density of the field 
stars. Integrating eq. ([3]) over the u-velocity leads to 
the density distribution of the induced polarization 
cloud. This has been calculated here without tak- 
ing into account th e self-gravity of the background 
medium. However, iFuchsl ( 2004 ) has shown that in 
linear approximation the effects of self-gravity are 
not important for the dynamics of the polarization 
cloud. 

3. POTENTIALS OF THE PERTURBING 
BODIES 

To start with, we consider in our analysis the po- 
tential of a point mass, which moves with the velocity 



vo along the y-axis 
*i = 



Gm 



Vx 2 + (y - v i) 2 + : 



(4) 



Its Fourier- Transform can b e calculated using for- 
mulae (3.754) and (6.561) of iGradshtevn fc Rhvzik 
(l2000h as 

= _^mJ_ e -& y v t 



2tt 2 fc 2 

Next, the potential (j4|) is generalized to 

Gm 



v/rjj + x 2 + (y - v i) 2 + : 



(5) 



(6) 



which corresponds to an extended body with the 
mass distribution of a Plummer sphere {p ex l/rg(l + 
r 2 /rg) - ^ ; Binnev fc Tremaine 1987t ) The Fourier 



transform of a moving Plummer sphere is given by 
Gm Tq 



2tt 2 k 



(7) 



where K\ denotes the modified Bessel function of 
the second kind. As third example we consider a 
perturber which has the mass density dist ribution of 
a Hern quist profile (p cx (ro/r)(ro + r)~ 3 : iHernquist 
Its gravitational potential is given by 



Gm 



ro 



(8) 



The Fourier transform of a moving Hernquist 
sphere can be calculated using eq. (3.722) of 
Gradshtevn fc Rhvzik! ^OOOl FI leading to 



$k = - 



Gm 1 
2^ fc 2 



1 + fcro cos(fcro)si(fcro) 



-kr sin(fcr )ci(fcr )]e- ikyVot , 



(9) 



where si and ci denote the sine- and cosine-integrals, 
respectively. The model of a Plummer sphere has of- 
ten been used in numerical simulations of the accre- 
tion and their eventual disruption of satellite galaxies 
in massive parent galaxies. Plummer spheres have 
constant density cores, whereas numerical simula- 
tions of the formation of galactic haloes in cold dark 
matter cosmology show t hat dark haloes may have 
a ce ntral density cusp ( Navarro. Frenk. fc White! 
19971 ). Thus models of a Plummer or a Hernquist 
sphere should encompass the range of plausible mod- 
els for satellite galaxies. The density in both models 
falls off radially steeper than found in the cold dark 
matter galaxy cosmogony simulations. This mimics 
the tidal truncation of satellite galaxies in the grav- 
itational field of their parent galaxies. 



3 We use the identity sinfcr : 
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cosfcr in eq. (3.722). 
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4. DYNAMICAL FRICTION 
The ensemble of stars is accelerated by the mov- 
ing perturber as 

< v >= - J d 3 x J d 3 v/(x, v)V$i , (10) 

where / denotes the full distribution function / = 
fa + fi- The contribution from /o cancels out, and 
introducing the Fourier transforms we find 



< v > = 



cfx / cfv / d 3 kik$ 



(11) 



From symmetry reasons the acceleration vector < 
v > is expected to be oriented along the y-axis. In 
cq. (jlip the frequency uj, and similarly uj' , is given 



according to eqns. ©,© an< l © by i 



-k y v - iX 



where we have introduced following Landau's rule a 
negative imaginary part, which we will let go to zero 
in the following. Moreover, = $_k so that the 
potential is a real quantity. Equation (jll[) simplifies 
to 



< v >= (2tt) 3 J d 3 v J d 3 kik$_ k / k e 2At . (12) 

Using expression (|3|) and taking the limit A — > we 
get 



< v >= 



(27r) 5 / 2 7rn b 



d 3 k|$ k | 2 ^k y v e~%^- . (13) 



The Fourier transform of any potential with spheri- 
cal symmetry depends only on k = |k|. Thus it fol- 
lows immediately from eq. (| 13[) that indeed the two 
acceleration components (< v x > = < v z > = 0) as 
anticipated. Only in the direction of motion of the 
perturber there is a net effect. According to New- 
ton's third law the drag force exerted on the per- 
turber is given by mv = -mi, < v y > where is 
the mass of a background particle, so that the drag 
force is anti-parallel to the velocity of the perturber. 
In order to evaluate the integrals over the wave num- 
bers in eq. (|13|) it is advantageous to switch from 
Cartesian form k x ,k y ,k z to a mixed representation 

k y , k = ^/k x + k y + k^ , arctan(k x /k z ), and we ob- 
tain for the deceleration the general result 
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In the case of a point mass formula ([5]) implies 
A = fc max /fc m j n . This result was first obtained 
in this form by Kama id ( 1972) and is identica l to 
Chandrasckhar's formula (jChandrasekharl 119431 ) . if 
m + mi, w m. The Coulomb logarithm diverges in 
the familiar way both on small and large scales, i.e. at 
k~l x and fc"^, respectively. The Coulomb logarithm 
of the dynamical friction force exerted on a Plummer 
sphere can be calculated by inserting eq. ([7]) into (fTS")) 
leading to 



InA 



dk kKf(r k). 



(16) 



If the Plummer radius rg shrinks to zero, expression 
(fTB)) changes smoothly into the Coulomb logarithm 
of a point mass, because lim ro _>o ro-Ki(feo&) = k^ 1 . 
The integral over the square of the Bessel functions 
in eq. (fT6|) can be eva l uated using formula (5.54) of 
Gradshtevn fe Rhvzikl (|2000l ). 



l nA = - JLjHlii [Kf(r k min ) - K (r k min )K 2 (r k min )] 

(17) 

which is approximately 



InA w —1/2 — ln(rofc m in) ; 



(18) 



in the limit of rofc m j n <^ 1. This modified Coulomb 
logarithm conve r ges o n small scales precisely as 
found by IWhitd (1970), but still diverges on large 
scales. A natural cut-off will be then the size of the 
stellar system under consideration. For a perturber 
with the density distribution of a Hernquist profile 
we find a Coulomb logarithm of the form 

r°° i 

InA = / dk — [1 + rofc cos(rofc)si(ro/c) 

— rofcsin(r fc)ci(r fc)] 2 . (19) 

It can be shown using the asymptotic expan- 
sions of the sine- and cos ine-integrals given by 



where erf denotes the usual error function. 
Coulomb logarithm is defined as 

4^4 rfc- 



InA = 



G 2 m 2 



dk fc^kl 



InA 

The 

(15) 



Abramowitz fc Stegunl (|1972r j that the integrand in 
expression (fT9|) falls off at large k as 4ro(rofc) -5 . 
Thus the Coulomb logarithm converges at small 
scales. This is expected because, although the den- 
sity distribution has an inner density cusp, the de- 
flecting mass 'seen' by a field star with a small im- 
(14)P ac "t parameter scales with square of the impact pa- 
rameter. At small wave numbers a Taylor expan- 
sion shows that the square bracket in expression (jT9j) 
approaches 1 so that we find a logarithmic diver- 
gence of the Coulomb logarithm as in the case of 
the Plummer sphere. In Fig. 1 we illustrate the 
Coulomb logarithms of the Plummer sphere and a 
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Fig. 1. Coulomb logarithms of the Phimmer sphere 
(solid line) , and a sphere with a Hernquist density profile 
(dashed line). Fhe dotted line indicates — ln(27rro/A ma x)- 
ro denotes the radial scale lengths of the spheres and 
Amax is the upper cut-off of the wavelength of the den- 
sity perturbations (see text). The inlet shows the cumu- 
lative mass distributions of the Plummer and Hernquist 
models. 



sphere with a Hernquist density distribution accord- 
ing to cqns. (fT7|) and (fT9j) as function of rofcmin- 
Fig. 1 shows clearly that at given mass small sized 
perturbers experience a stronger dynamical friction 
force than larger ones. Since the cut-off of the wave 
number at fc m in is determined by the radial extent 
of the stellar system, which corresponds roughly 
to one half of the largest subtended wave length 
Amax = 27r/fc min , we use actually ro/(A roax /4) as ab- 
scissa in Fig. 1. For comparison we have also drawn 
— ln(rok m i n ) in Fig. 1. The insert shows the cumu- 
lative mass distributions of both mass models. The 
half mass radius of the Hernquist model measured 
in units of ro is about twice as that of the Plummer 
sphere. Thus for a proper comparison of the drag 
forces exerted on a Plummer sphere and a sphere 
with a Hernquist density profile the dashed line in 
Fig. 1 should be stretched by a factor of about 2 to- 
wards the right. But it is clear from Fig. 1 that the 
drag force exerted on a Plummer is always larger that 
the drag on a sphere with a Hernquist density profile. 
This is to be expected because of its shallower den- 
sity profile. The logarithm — ln(rok m i n ), although 
being the asymptotic expansion of the Coulomb log- 
arithms (fTT)) and (flU)) for r k min — > 0, is not a good 
approximation at larger rpfc m i n . There is a system- 



atic off-set relative to the true Coulomb logarithms 
which is given explicitely in eq. (fl8|) for the case of 
the Plummer sphere. It could be argued that per- 
turbers can be deformed by tidal fields. However, its 
effect is expected to be small (a higher order effect). 
There is a further effect if the perturbers are gravita- 
tionally bound systems themselves like globular clus- 
ters or dwarf satellite galaxies. Such objects can and 
do loose mass due to tidal shocking and both effects 
can even be of comparable magnitude. Finally, we 
want to mention that our analysis can be extended, 
in a straightforward way, to anisotropic velocity dis- 
tributions of the field stars. iFuchs fc Athanassoula 



(|200a ) have shown that the velocity dispersion in the 
solution of the Boltzmann equation ([3]) is replaced 
by an effective velocity dispersion which depends on 
the semi-axes of the velocity ellipsoid and its orien- 
tation relative to the wave vector k. We intend to 
present our results in an forthcoming paper (Esquivel 
& Fuchs, in preparation). 

5. ACKNOWLEDGEMENTS 

O.E. gratefully acknowledges financial support 
by the International-Max-Planck-Research-School 
for Astronomy and Cosmic Physics at the Univer- 
sity of Heidelberg. ] 

REFERENCES 

Abramowitz M., Stegun LA. 1972, Handbook of Mathe- 
matical Functions, Dover, New York 

Binney J., Tremaine S., 1987, Galactic Dynamics, Prince- 
ton University Press, Princeton 

Chandrasekhar S. 1943, ApJ, 97, 255 

Dekker E., 1976, Phys. Rev. D, 24, 315 

Fuchs, B. 2004, A&A, 419, 941 

Fuchs B., & Athanassoula, E. 2005, A&A, 444, 455 

Gradshteyn, I. S., & Rhyzik, I.M. 2000, Table of In- 
tegrals, Series, and Products, Academic Press, New 
York 

Henon, M. 1973, in: Dynamical structure and evolu- 
tion of stellar systems, G. Contopoulos, M. Henon, 
D. Lynden-Bell (eds.), Lectures of the 3rd Advanced 
Course of the Swiss Society of Astronomy and Astro- 
physics, Geneva Observatory, Sauverny, p. 182 

Hernquist, L. 1990, ApJ, 356, 359 

Kalnajs, A. J. 1972, in: Gravitational N-Body Problem, 
M. Lecar (ed.) Reidel Publ. Comp., Dordrecht, p. 13 
Lynden-Bell, D., & Kalnajs, A. J. 1972, MNRAS, 157, 1 
Marochnik, L. S. 1968, Soviet Ast., 11, 873 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, 

ApJ, 490, 493 
Stix, T. H. 1992, Waves in Plasmas, AIP, New York 
Tremaine, S., & Weinberg, M. D. 1984, MNRAS, 209, 
729 

White, S. D. M. 1976, MNRAS, 174, 19 



